Canonical form of master equations and characterization of non-Markovianity 
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Master equations govern the time evolution of a quantum system interacting with an environment, 
and may be written in a variety of forms. Markovian master equations, in particular, can be cast 
in the well-known Lindblad form. Any time-local master equation, Markovian or non-Markovian, 
may in fact also be written in Lindblad-like form. A diagonalisation procedure results in a unique, 
and in this sense canonical, representation of the equation. This representation may be used to 
fully characterize the non-Markovianity of the time evolution. Recently, several different measures 
of non-Markovianity have been presented. Their common underlying definition of non-Markovianity 
is whether negative decoherence rates may appear in the Lindblad-like form of the master equation. 
We therefore propose to use the negative decoherence rates themselves, as they appear in the unique 
canonical form of the master equation, as a primary measure to more completely characterize non- 
Markovianity. The advantages of this are especially apparent when many decoherence channels are 
present. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc 



An open quantum system is a quantum system whose 
dynamics is determined both by interactions internal to 
the system, and by influences from an environment. As 
no physical system is truly isolated, this is a situation 
which applies very widely. Markovian or memory-less 
behavior leads to master equations in the so-called Lind- 
blad form 0, Lindblad master equations have been 
extensively used to describe phenomena in e.g. quantum 
optics, semiconductor physics and atomic physics, rang- 
ing from the decay of an atom to a quantum- mechanical 
description of Brownian motion [H, H| . 

The non-Markovian case is less well understood, but 
is becoming increasingly relevant as our ability to ex- 
perimentally control quantum systems develops. Sim- 
ple examples are a damped harmonic oscillator and a 
damped driven two-level atom Q. Also, the vast major- 
ity of calculations of error thresholds for quantum com- 
puting makes the assumption that the noise processes are 
Markovian. This is not necessarily physically realistic Q . 

Recently, several different measures of non-Markov- 
ianity have been proposed They all have in com- 
mon that they take as the underlying definition of (possi- 
bly time-dependent) Markovianity that the master equa- 
tion may be written in a local-in-time Lindblad-like form, 
to be defined below, where all decoherence rates jk (t) are 
positive at all times. Non-Markovian time evolution is 
characterized by the occurrence of negative decoherence 
rates in the master equation. The measure in Q is based 
on whether it is possible for the trace distance between 
two initial states to increase as a function of time. If the 
time evolution starting from a time to is completely posi- 
tive (CP), then the trace distance cannot exceed its inital 
value at to, but it could first decrease and then increase 
again. This would be a signature of non-Markovian be- 
havior. The measure in is similarly based on whether 
entanglement with an environment can increase, Q con- 



siders the flow of Fisher information, and [fj the amount 
of isotropic noise that has to be added to make the time 
evolution (time-dependent) Markovian. 

In order to obtain a fundamental characterisation un- 
derlying the many different measures proposed, it would 
seem natural to focus on their common point of origin. A 
primary measure of non-Markovianity would then be the 
negative rates 7&(i) themselves. Alternatively, one could 
recognise that any function of the 7fc(i) that captures 
their negativity in an adequate way (exactly how might 
vary for different contexts) will be a valid measure of 
non-Markovianity. The measures in [6 9] can be viewed 
in this way, and one may suggest other such functions. 
Individual measures such as increasing trace distance or 
increasing entanglement are still useful e.g. as experi- 
mentally accessible indicators of the presence of negative 
decoherence rates. If only one decoherence channel is 
present, they can be shown to be equivalent to knowl- 
edge of 7(£) @, [3 , but in general they neither require 
nor give full knowledge of the master equation. An ad- 
vantage with using the negative jk(t) themselves is that 
one obtains a unique measure that gives a more complete 
picture of non-Markovianity even when there are several 
decoherence channels present. Also, we do not need to 
solve the master equation to obtain the time evolution of 
the system, or optimise over initial states. 

If one is to base a measure of non-Markovianity di- 
rectly on the negativity of decoherence rates in a master 
equation, it is important that one considers a unique form 
the equation, since each master equation may be written 
in many ways. This is indeed a reason cited in Q for 
not basing a definition of non-Markovanity directly on a 
master equation. Fortunately, there is such a unique and 
canonical form. This result is a straightforward extension 
of the treatment in . A general form of time- local mas- 
ter equations has been used previously e.g. in [lol [llj ]. 
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but the existence of a unique diagonal form of the equa- 
tion has as far as we are aware not been emphasized. 
Neither is the generality of time-local master equations 
as widely recognised as it should be. We will therefore 
start by presenting a derivation of the canonical form of 
a time-local master equation, before discussing how this 
gives rise to natural ways to quantify non-Markovianity. 

I. FORMS OF MASTER EQUATIONS 

Under fairly general conditions, a master equation for 
p takes the form fl2l . [l3j 

p{t) = - l -[H,p(t)] + j\ s Ms)]ds. (1) 

Here H is the system Hamiltonian, and the memory ker- 
nel K, St t is a linear map describing the effects of the en- 
vironment on the system. The Born-Markov approxima- 
tion amounts to approximating the memory kernel by 
fcs,t[p(s)] ~ JCS(t — s)p(s). In this case, there is no mem- 
ory, and the system dynamics is Markovian. This leads 
to a master equation in Lindblad form [l], Q , 

p-^n-E| { L i L *p + p L i l » - 2L kpL[) , (2) 

where jk are positive decoherence rates, the operators 
Lk describe the different decoherence channels, and the 
unitary part of the time evolution corresponding to the 
Hamiltonian term may include effects arising from the 
environment. 



A. Time-local master equations 

Even for the case of non-Markovian dynamics, when 
the system 'remembers' its previous evolution, it has been 
shown via the time-convolutionlcss projection operator 
method [HI, EH that a general master equation ([T]) often 
can be written in a local-in-time form, p(t) = A t [p(t)], 
where A t is, among other requirements, a linear map 
such that At(p) is Hcrmitian and traceless for all p. For 
non-Markovian evolution, A f is timc-dcpcndcnt, and it 
is difficult to characterize what forms it may take in or- 
der to be compatible with physical or completely positive 
evolution. A description in terms of time-local master 
equations is not only convenient, but also seems to be 
necessary for constructing a quantum trajectory unrav- 
elling in terms of possible 'paths' the time evolution of 
the system may take [TJ I16H19} ■ 

As an alternative simple proof [2(| that memory-kernel 
master equations can in fact be written in time-local 
form, consider some evolution process described by the 
linear (usually completely positive) map p(t) = <fi t [p(0)], 
which satisfies a memory-kernel master equation of the 
general form ([T]) above. If it is assumed that the map 



4> t is invertible for the time interval considered, i.e., that 
there exists a linear map <j)^ 1 satisfying (j)^ 1 o <p t = I 
where / is the identity map, then one can write 

p(t) = f ds(K a , t o(f> a )\p(p)] 
Jo 

= A t [p(t)], (3) 

where A t := J * ds K. s j o <p s o tj)^ 1 , which is of time-local 
(albeit in general complicated) form, and the Hamilto- 
nian evolution is absorbed into K.. 

The assumption that the evolution is invertible is vi- 
olated only if two initially distinct states evolve to the 
same state at some time t (e.g., if an 'equilibrium state' 
is reached within a finite time). As shown in [20j . even 
in this case it is sometimes possible to describe the evo- 
lution by a time-local master equation. The restriction 
that the time evolution be invertible need not be serious. 
Physically relevant examples where the time evolution is 
invertible for any finite time include spontaneous emis- 
sion and optical phase diffusion [2l[ . The latter is also an 
example of a time-local master equation where the time 
evolution is not of standard Lindblad form. Also, even if 
the time evolution is not invertible for isolated points in 
time, but the master equation exists for other times, one 
can still use it for a characterization of non-Markovianity. 
Typically, decoherence rates would approach infinity for 
times when the evolution is not invertible (20l . I22I [23j . 

B. Lindblad- type canonical form of master 
equations 

We will now show that any time-local master equa- 
tion can be cast in Lindblad-like form, where some of 
the time-dependent decoherence rates jk (t) may become 
negative for some time intervals. The derivation is a 
straightforward generalisation of the results in Q , but de- 
serves to be more generally known. Lindblad-like terms 
with the 'wrong sign' are related to interaction between 
the environment and the system in such a way that 
the system may recohere, reversing earlier decay pro- 
cesses mmm. 

A general local-in-time master equation has the form 
p = J2Mt)pBl(t). (4) 

From the requirements that p remain Hermitean for all 
time and that the trace of p be preserved, we can ex- 
tract some general properties regarding the form of the 
master equation. For a state space of dimension d, we 
first introduce a complete set of N(= d 2 ) basis operators 
{G m ; m = 0, 1, 2, . . . N} with the properties 

G = l/Vd; G n = Gl; Tx[G m G n ] = 6 mn , (5) 
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where / is the identity. The last condition implies, by 
putting m = 0, that Tr[G„] = for n ^ 0. 

For brevity, we will suppress the time dependence in 
quantities below, but everything except the basis opera- 
tors G n may be time dependent. We can expand 



We now observe a crucial feature of this last result that 
seems not to have been appreciated before, but which en- 
ables us to derive a main result of this paper. Thus, we 
take advantage of the Hermitean nature of the decoher- 
ence matrix to write it in diagonal form, 



A k = G ia%k , B k = ^ Gj bji 



dij = ^ U ik j k U, 



so that p = J2ijJ2 k a ikb*j k GipGj. If we now define the 
quantities aj = J2k a ikb*j k , then 

P = ^2 CijGipGj. 

ij 

Using the fact that p and hence p are Hermitean, then 
^ ' CijGipGj = 'S ] c*jGjpGi = 'S ' CjiGjpGi, 



so that Cij = c^. Thus the Cy are the elements of an 
N x N Hermitean matrix. By separating out the i = 
and j = terms, we can write the expression for p as 



where the Ut k are the eigenvectors of d, J2 k Ui k U* k = 5ij, 
and the eigenvalues ^ k of d are real but not necessarily 
positive at all times. Define the time dependent operators 

Lfc(t) = J2U ik {t)Gi. 
We then have, via ([7]), the canonical form 

JV-l 

p = -i[H(t),p] + Y, 7fc(*) \L k {t) P L\{t) (8) 



\L\{t)L k {t)p-\pLl{t)L k {t) 



Vd I \ \fd 



This may be recognised as being of Lindblad-type 
= ~^fP Jr \ Z-^i 7/7i < ^ i ) | ~Jl*~*i I dijGipGj, f orm except that the tracclcss operators L k (t) are 

time dependent in general, as are the eigenvalues j k (t) 
which also arc not guaranteed to be positive. The 
decoherence channels arc orthogonal in the sense that 
Tr(L^L n ) = S mn . More general Lindblad-type forms 
correspond to non-diagonal decompositions of d. The 
consequences for ^ k (t) and L k (t) of the requirement that 
the time evolution should be completely positive, or even 
positive, have not been considered here. Thus, for in- 
stance, there is no guarantee that an arbitrary master 
equation of this form will yield only positive eigenvalues 
for p for all time. Note that under any unitary transfor- 
mation p — > V^pV^Y , e.g. to an 'interaction' picture, 
both H{t) and L k (t) will change. However, the j k (t) are 
invariant under any such change of picture f25j . 



where now the quantities dij(= cy) will be elements of 
an (N — 1) x (N — 1) 'decoherence' Hermitean matrix d. 
Furthermore, writing 

n _ 1 coo v-^ c l0 „ 
and since Cqo is real and cqj = c* , we have 



p = Cp + pC f + Y dijGipGj. 



(6) 



If we now apply the condition that Tr[p] = we get 



Tr 



which tells us that C + = — Y^,ijM '!,,(■' ,(■',. The mas- 
ter equation is therefore, with H = —2i(C — C') t 



P=2 [{C-C^)p + p(C^ -C) 



(7) 



-(C + C^)p + p(C^ + C)] + J2 d ij G iP G j 



[H,p] + ]T dtf(i) (g x pG, - l -G 3 G lP - l -pG^, 



This is the kind of structure obtained in 0, though as 
they are considering quantum semigroups, their decoher- 
ence matrix d is independent of time, as is H. 



II. CHARACTERIZING NON-MARKOVIANITY 

Time-local master equations are thus widely applica- 
ble, and any time-local master equation may be rewritten 
in the canonical form in ([5]) . The at most N — 1 decoher- 
ence rates 7fc(t) and operators L k (t) in ((SJ) arc obtained 
from a diagonalisation procedure and consequently the 
7fc(i) are unique, with the L k (t) unique up to degenera- 
cies and choice of basis. If the j k (t) are positive at all 
times, then the time evolution is completely positive in 
any time interval. This may be termed 'time-dependent 
Markovian' evolution, as done e.g. in The defining 

feature of non-Markovianity is then whether any of the 
7fc(i) may become negative. It is therefore natural to use 
the functions 



/*(*) 



(9) 
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to characterize the non-Markovianity of the time evo- 
lution. Each fk(t) describes the non-Markovanity in 
an individual decoherence channel. We can either use 
fk(t) as they are or define a function of them, such as 

their sum, f(t) = J2k=i /fe(*)i ^ we wish, normalised 
by some function of N. Another possibility is to use 
Fk(t) = Jq fk(t)dt, to characterize the 'total amount of 
non-Markovianity in channel k up until the time £', or 
Fk(t) for the 'total amount of non-Markovianity 
up until the time t\ Either of these quantities may be 
normalized e.g. by the total elapsed time t. 

We may also define a discrete measure of non- 
Markovianity as the number of negative 7fc(t), a 'non- 
Markov index'. In the sense that Tr(Lj„L„) = 6 mn , 
the non-Markovian part of the dynamics takes part in 
a region of 'evolution space' which is orthogonal to the 
Markovian region. For example, for a two-level system, 
if Markovian behavior is taking place along the x direc- 
tion of the Bloch vector, then non-Markovian behavior 
can take place along the y and z directions. The non- 
Markov index characterizes the dimension of the space 
of non-Markovian evolution. 

To illustrate the importance of using the canoni- 
cal form of the master equation to characterize non- 
Markovianity, consider the Lindblad-type equation 

p = [2j{t) + j(t)] \2a x pa x + 2a y pa y - Ap] 
— j(t) [2(j-pa + — o + <j-p — p<j + oJ\ 
—j(t) [2a+p<7- — o-<j + p — /5<7_<7_|_] . (10) 

This may at first sight look non-Markovian for j(t) > 
or 27(i) + j(t) < 0, but the equation is not in canonical 
form and may in fact be rewritten 

P = [l{t) + 7(t)] [2a x pa x + 2a y pa y - Ap] , (11) 

which is Markovian for j(t) + j(t) > 0. 

In the simplest case with only one non-zero 7fc(i) = 
7(t) and time independent Lk, the trace distance 
measure of non-Markovianity becomes equivalent to 
— j(t) exp[— L dt'^it')] and the entanglement-based 
measure becomes equivalent to —2 J u\ <0 j{t)dt [8|. The 



isotropic noise one must add in order to make the time 
evolution completely positive is also directly implied by 
min[7(i), 0)]. At times when the time evolution is not 
invertible, some 7fc(t) in the time- local master equation 
become singular. So do the trace-distance and entangle- 
ment measures of non-Markovianity. 

When there is more than one non-zero 7fc(i), that is, 
more than one decoherence channel is displaying non- 
Markovian behavior, a single number cannot capture the 
full nature of the non-Markovianity. One may consider 
e.g. the trace distance between not just one pair of states, 
but instead the time-dependent pairwise trace distances 
between a selection of initial states. If the underlying 
definition of non-Markovianity is the behavior of the neg- 
ative rates 7fc(t) themelves, however, then any 'complete' 
set of measures for non-Markovianity must by definition 
be sufficient to reconstruct the 7fc(i). It may be that 
such a complete set of e.g. trace distances, together with 
knowledge of the respective initial states, is sufficient to 
reconstruct not just the jk(t) but also the Lk(t). and 
thus the complete master equation. The relation between 
7fc(i) and the measures in @-@| for time dependent H(t) 
and Lk(t) remains to be investigated. The invariance of 
7fc(i) for different Hk{t) and Lk(t) may be an advantage. 

For a full description of non-Markovianity, one should 
include also information about H(t) and Lk(t). Although 
7fc(t) are unique and invariant under basis transforms, 
different H(t) and/or Lk(t) for the same 7fc(£) will corre- 
spond to different relaxation behavior. One might base 
a full description of non-Markovianity on the 'interac- 
tion picture' where H{t) = 0. Complete positivity will 
similarly also depend on the full dynamics, when the evo- 
lution is non-Markovian. 

In summary, there exists a canonical time-local form 
for master equations, valid both for Markovian and non- 
Markovian time evolution. If non-Markovian time evolu- 
tion corresponds to negative time dependent decoherence 
rates, then the decoherence rates in the canonical form 
of a master equation may be used to describe the non- 
Markovian character of the time evolution. 
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